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The structure and function of the plant microbiome is driven by plant species and 
prevailing environmental conditions. Effectuated by breeding efforts, modern crops 
diverge genetically and phenotypically from their wild relatives but little is known 
about consequences for the associated microbiota. Therefore, we studied bacterial 
rhizosphere communities associated with the wild beet B. vulgaris ssp. maritima grown 
in their natural habitat soil from coastal drift lines (CS) and modern sugar beets 
(Beta vulgaris ssp. vulgaris) cultivated in CS and potting soil (PS) under greenhouse 
conditions. Analysis of 16S rRNA gene fingerprints and pyrosequencing-based amplicon 
libraries revealed plant genotype- and soil-specific microbiomes. Wild beet plants harbor 
distinct operational taxonomic units (OTUs) and a more diverse bacterial community 
than the domesticated sugar beet plants. Although the rhizospheres of both plant 
genotypes were dominated by Proteobacteria and Planctomycetes, 37.5% of dominant 
OTUs were additionally detected in the wild beet rhizosphere. Analysis of the cultivable 
fraction confirmed these plant genotype-specific differences at functional level. The 
proportion of isolates displayed in vitro activity against phytopathogens was lower for 
wild beet (<45.8%) than for sugar beet (<57.5%). Conversely, active isolates from the 
wild beet exhibited stronger ability to cope with abiotic stresses. From all samples, 
active isolates of Stenotrophomonas rhizophila were frequently identified. In addition, 
soil type-specific impacts on the composition of bacterial communities were found: 
Acidobacteria, Chloroflexi, and Planctomycetes were only detected in plants cultivated 
in CS; whereas Bacteroidetes and Proteobacteria dominated in PS. Overall, in comparison 
to modern sugar beets, wild beets were associated with taxonomically and functionally 
distinct microbiomes. 
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INTRODUCTION 

The domestication of plants began with the identification of wild 
plant species exploitable for food, animal feed, or other domestic 
purposes. Selective crop breeding programs are mainly oriented 
toward increasing food production and high-yielding varieties, 
percentage of usable plant parts, and resistance against diseases 
(Akhalkatsi et al, 2012). The genetic background of sugar beet 
is narrower compared to other crops because the sugar beet 
breeding had started in the late eighteenth century when lines 
accumulating sugar in the storage root were selected from crosses 
made with chard and fodder beet (Fischer, 1989). In contrast to 
many other crops, the ancestor of domesticated sugar beets (Beta 
vulgaris ssp. vulgaris L.) is known and still distributed throughout 
Europe, especially in the Mediterranean region along coastal drift 
lines. The genome of the wild form B. vulgaris L. ssp. maritima 
(L.) Arcang. (trivial name: wild or sea beet) differs significantly 
from domesticated sugar beet lines (Dohm et al., 2014). For 



example, it shows high variation at the vernalization 5-locus, 
and therefore most of the plants are perennials in contrast to all 
cultivated lines, which have a biennial life cycle enabling seed pro- 
duction (Dohm et al., 2014). Today, modern sugar beet cultivars 
that accumulate enormous amount of sugar are of importance, 
not only for commercial sugar production, but also as renew- 
able resource, e.g., for bioethanol production (Dodic et al, 2009). 
Although several studies analyzing sugar beet-associated bacterial 
communities already exist (Zachow et al, 2008; Mendes et al., 
2011; Shi et al, 2014), little is known about the relation of the 
microbiome to existing cultivars and their wild relatives. 

Plants have recently been recognized as meta-organisms, 
due to a close symbiotic relationship with their microbiome. 
Comparable to humans and other eukaryotic hosts, plants also 
harbor a "second genome" that fulfills important host functions 
including protection against biotic and abiotic stress (Berendsen 
et al., 2012; Berg et al., 2013). Interestingly, plant-associated 
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bacteria derived from different origins. Microorganisms can be 
transmitted by pollen and seeds (Ftirnkranz et al., 2012); the lat- 
ter was also shown for sugar beet (Dent et al., 2004). However, 
specific microorganisms are enriched from the surrounding envi- 
ronment as well, e.g., attracted by root exudates containing 
carbohydrates, proteins, and vitamins (Chaparro et al., 2013). 
In result of these processes, each plant harbors to a certain 
degree specific microbes (rev. in Berg and Smalla, 2009; Bulgarelli 
et al, 2012). This specificity was also shown for the plant- 
associated microbiome at cultivar level, e.g., for rice (Engelhard 
et al., 2000; Hardoim et al, 2011) as well as for maize (Haichar 
et al, 2008; Philippot et al., 2013). The evolutionary relationship 
in wheat-microorganism interactions was revealed already by 
Germida and Siciliano (2001); ancient wheat cultivars were colo- 
nized by phylogenetically diverse rhizobacterial isolates, whereas 
the rhizosphere of modern cultivars was dominated by fast- 
growing Proteobacteria. Therefore, we developed the hypothesis 
that the microbiome of the wild beet harbors a high degree of 
specific microorganisms in comparison to the sugar beet crop. 
We expected this degree of specificity also for the functions of 
the associated bacteria. The wild beet plants commonly grow in 
salinated dune soil of coastal drift lines. Thus, protection against 
abiotic stress by the microbiome should be more important than 
against biotic stress because plant diseases are not known for these 
populations (Biancardi et al., 2012). As reason for the specific 
microbiome of the wild beet, we assumed not only the genetic 
background of the wild ancestor; we also hypothesized an impact 
of the microbiome of the coastal drift line soil. 

The objective of this study was: to determine key players in the 
rhizosphere microbiomes of (i) wild beet plants grown in coastal 
drift line soil (WB-CS), and (ii) domesticated sugar beets cv. 
BERETTA (SB) in dependence of the soil type (coastal drift line 
soil, CS, and potting soil, PS), and (iii) to determine the poten- 
tial of isolates to cope with biotic and abiotic stresses in successive 
screenings. 

MATERIALS AND METHODS 
SAMPLING AND EXPERIMENTAL DESIGN 

Wild beet plants (WB) with the synonymous common name 
sea beet [B. vulgaris ssp. maritima (L.) Arcang.] were sampled 
from the drift line at the Mediterranean Sea coast in Slovenia (N 
45.590725, E 13.719469). Along a 200 m coastal drift line, four 
independent samples were taken. On average, plants were 20 cm 
high. Leaves were unfolded and separately spread on the beet 
heads without flower buds or seed development. Additionally, we 
collected coastal drift line soil (CS, pH 9.5, electric conductance 
S = 1728 u,S cm -1 ) from the layer of 0-10 cm depth. Samples 
were placed into sterile plastic bags and transported to the lab- 
oratory at 4°C. Commercial sugar beet (SB) seeds cv. BERETTA 
KWS were provided by KWS SAAT AG (Einbeck, Germany). 
Seeds (10 per pot) were planted in a potting bulk soil (PS, 
soil:sand:vermiculite mixture 3:1:1, potting soil basis: "Gramofior 
Profi-Substrat-Topfpikier M+Ton+Fe," GBC, Kalsdorf, Austria, 
pH (CaCh) 5.8, electric conductance S = 1330 |xS cm -1 , con- 
taining white and black peat with 90 kg m~ 3 moist clay, 1.0 kg 
m~ 3 PG-Mix (Greenworld, Wels, Austria), 50 g m~ 3 Radigen® 
(Terrafior, Iserlohn, Germany); sand: "Maxs Spielsand®," Scherf 



GmbH & Co KG, Hartberg, Austria; vermiculite: "Verm 3-6 mm," 
Ratioform, Vienna, Austria) as well as in coastal bulk soil from 
the drift line (CS). For SB-PS, five replicates, and for CS — due 
to the limited amount of available soil and plant material — 
three (SB-CS) replicates were incubated under greenhouse con- 
ditions (12 h light/dark at 18°C, 60% humidity) for 2 weeks. 
To isolate microorganisms, 20 g of roots with adhering soil were 
mixed with 50 ml sterile 0.85% NaCl for 3 min in a laboratory 
blender (BagMixer, Interscience, Mourjou, France). Suspensions 
were used for cultivation-independent and -dependent analy- 
ses. Samples were labeled as follows: CS — coastal drift line soil, 
PS — potting soil, SB-CS — sugar beet plants cultivated in coastal 
drift line soil, SB-PS — sugar beet plants cultivated in potting soil, 
WB-CS — and wild beet plants grown in coastal drift line soil. 

MOLECULAR FINGERPRINTS OF 16S rRNA GENES 

For cultivation-independent analyses, 8 ml of the above men- 
tioned suspension were centrifuged at high speed (16,000 x g, 
4°C) for 20 min and the resulting microbial pellets were stored 
at — 70° C until further processing. Total DNA of the rhizo- 
sphere communities were extracted by mechanical disruption and 
homogenization of the pellet using FastDNA Spin Kit for Soil and 
a FastPrep Instrument (MP Biomedicals, Illkirch, France) for 30 s 
at 5.0 m s _1 . DNA was additionally purified by the GeneClean 
Turbo Kit (MP Biomedicals, Illkirch, France) containing guani- 
dine thiocyanate to remove humic acids. Extracted DNA was 
treated with RNase (0.02 ng |xl _1 ) for 5 min at 65° C to obtain 
the template for PCR amplification of 16S rRNA genes from total 
community DNA. 

Single strand conformation polymorphism (SSCP) analysis 
was carried out according to Schwieger and Tebbe (1998) using 
a 8.0% polyacrylamide gel running 26 h for bacterial com- 
munity analysis at 400 V. SSCP-PCR was performed accord- 
ing to Koberl et al. (2011). Gels were transmissively scanned 
(Epson perfection 4990 Photo, Nagano, Japan) to obtain dig- 
itized gel images. Normalization and cluster analysis of band 
patterns, evaluated on band intensity, were carried out with the 
GelCompar II program (Applied Maths, Sint-Martens-Latem, 
Belgium). The cluster analysis was performed using the follow- 
ing settings: dendrogram type: unweighted pair group method 
with arithmetic mean (UPGMA); curve based similarity coeffi- 
cient: Pearson correlation; optimization 4%, and position tol- 
erances 1%. Background correction was applied for each track. 
The Pearson correlation index for each pair of lanes within a gel 
was calculated as a measure of similarity between the community 
fingerprints. 

AMPLIC0N PYROSEQUENCING AND BI0INF0RMATIC ANALYSIS 

The hypervariable V4-V5 region of the 16S rRNA gene 
(Escherichia coli positions 515-927) was amplified in a nested 
PCR approach for 454 pyrosequencing to analyze the taxonomic 
composition of the bacterial rhizosphere community. The first 
PCR was conducted with the primer pair 27f/1492r (Lane, 1991), 
while the second PCR targeted the V4-V5 region with the primer 
set that contained the 454 pyrosequencing adaptors and sample- 
specific tags (Table SI). The reaction mixture for the first PCR 
(10|xl) contained 1 x Taq&Go (MP Biomedicals, Eschwege, 
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Germany), 0.1 |xM of each primer and 1 (il of template DNA 
(95°C for 5min, 30 cycles of 95°C for 30 s, 57°C for 30 s, 72°C 
for 90s, and elongation at 72° C for 5min). The second PCR 
(60 (il) was performed by using 1 x Taq&Go, 1.5 mM MgCi2, 
0.1 [iM of each primer and 3 |xl of first PCR template (95° C 
for 5min, 32 cycles of 95°C for 20 s, 54°C for 15 s, 72°C for 
30 s, and elongation at 72° C for lOmin). PCR products of four 
replicates per samples of the same habitat were pooled in equal 
molarities and purified with the Wizard™ SV Gel and PCR Clean- 
up System (Promega, Madison, USA). 16S rRNA gene ampli- 
cons were pyrosequenced with Roche 454 FLX GS conducted 
by MWG Biotech (Ebersberg, Germany). De-multiplexed raw 
sequences were processed using the open source software package 
QIIME (release 1.8.0; Caporaso et al., 2010). Prior to denoising 
(Reeder and Knight, 2010) sequences were quality (minimum 
average quality score in reads: 25) and length filtered (430- 
450 bp). Chimeric sequences were detected via ChimeraSlayer 
and subsequently removed. Chimera check was followed by 
excluding plant-originated plastid sequences using BLASTn algo- 
rithm (Altschul et al., 1990). Remaining sequences were clus- 
tered at 97% similarity using the UCLUST algorithm (Edgar, 
2010) and taxonomic assignment of representative sequences 
were performed using the RDP naive Bayesian rRNA classifier 
(Wang et al., 2007) based on the reference database Greengenes 
release 13_5 (De Santis et al., 2006). Read number per sam- 
ple was normalized to 5578. Rarefaction analysis (species level 
at 97%), richness estimates, and diversity indices were calcu- 
lated; Shannon (1997) and Chaol (Chao and Bunge, 2002) 
indices were calculated based on the complete linkage cluster- 
ing data. For OTU-based analysis, only OTUs accounting for at 
least 1% of total reads were considered. Classifications of the 
reads (one representative per OTU) were performed by man- 
ual alignment of representative sequences with 16S rRNA refer- 
ence gene sequences from NCBI database using BLASTn algo- 
rithm. Raw pyrosequencing data were deposited at the National 
Center for Biotechnology Information under the BioProject num- 
ber PRJNA233435 with the SRA accession numbers SRX652486 
(WB-CS), SRX652836 (CS), SRX652838 (SB-CS), SRX652839 
(SB-PS), SRX652840 (PS). 

ISOLATION AND CHARACTERIZATION OF BACTERIA 

The homogenized suspensions of wild beet rhizosphere repli- 
cates were used for dilution and plating on R2A (Roth, Karlsruhe, 
Germany), and Kings B amended with ampicillin (50 (ig ml -1 ), 
novobiocin (45 (ig ml -1 ), and cycloheximide (50 |xg ml -1 ). Each 
dilution was plated in duplicates. Plates were incubated for 5 
days at room temperature (RT) and colony forming units (CFU) 
were counted to calculate the means of colonies (login CFU) 
based on fresh weight (fw). To obtain the microbial communi- 
ties of commercial sugar beet plants (SB-CS, SB-PS), 2 weeks 
after germination the culturable fraction of the rhizospheres 
were harvested, suspended, diluted, and plated as described 
above. If possible, 24 bacterial isolates were selected randomly 
for each replicate and subsequently cultured on nutrient agar 
(NA). Due to the limited amount of CS, fewer replicates were 
available and therefore, a lower number of bacteria were iso- 
lated. The isolates were purified and then stored at — 70° C in 



nutrient broth II (NB II) (Sifin, Berlin, Germany) stocks con- 
taining 15% (v/v) glycerol. Isolates were encoded using a com- 
bination of letters and numbers indicating: (1) plant/soil type 
(WB — wild beet plants grown in coastal drift line soil WB-CS, 
VN — sugar beet plants cultivated in coastal drift line soil, VS — 
sugar beet plants cultivated in potting soil), (2) replicate (1-4), 
(3) consecutive number of the isolate per replicate, and (4) ori- 
gin of the medium (no further indication — bacteria from R2A, 
Ps — Pseudomonadaceae from Kings B). Bacteria with the high- 
est antagonistic activity and stress tolerance were identified by 
sequencing partially the 16S rRNA gene using the primer pair 
27f/1492r (Lane, 1991). Purified fragments were sequenced using 
LGC Genomics GmbH sequencing service (Berlin, Germany) and 
identified as described. Sequences obtained were submitted to 
EMBL Nucleotide Sequence Database under accession numbers 
KI024636 to KJ024701. Co-occurrence of 16S rRNA genes from 
isolates and OTUs obtained by pyrosequencing was determined 
by aligning corresponding regions using BLASTn algorithm. 

SCREENING FOR BACTERIA ANTAGONISTIC TOWARD PLANT 
PATHOGENS 

The antagonistic potential of randomly selected isolates was 
assessed by in vitro inhibition of sugar beet-pathogenic fungi 
Alternaria alternata Nees, Botrytis cinerea Persoon, Rhizoctonia 
solani Kiihn AG2-2IIIB, Sclerotinia sclerotiorum Fuckel, and 
Verticillium dahliae Klebahn V25 according to Berg et al. (2002). 
Zones of inhibition were measured 3-7 days after incubation 
at25°C. 

SCREENING FOR STRESS TOLERATING BACTERIA 

Overnight cultures grown in 10 ml NBII were used as inoculum 
(5 (il) for all bacterial assays. For desiccation assays, 20 |xl of cul- 
ture fluids were dried under sterile conditions. After 3, 6, 9, 16, 21, 
and 56 days, bacteria were re-suspended with the same volume 
NBII and dropped on NA. Re-cultivated cells were registered in 
a positive/negative response. For osmolarity stress, bacteria were 
cultivated in a 96-well plate filled with 145 |xl per well modi- 
fied Luria Bertani (per liter: peptone 3 g, meat extract 5 g) with 
various sodium chloride concentrations (1% steps, 0-20%) and 
incubated at 30° C for 48 h under agitation. In reactive oxygen 
species stress assays, bacteria were cultivated in NBII amended 
with different tellurite (1, 3, 5, 7, and 9mg ml -1 ) and H2O2 
(0.1, 0.3, 0.5, 0.7, and 0.9 mM) concentrations. The assays were 
performed in 96-well-plates filled with 195 (il medium per well 
and 5 (il 1:100 diluted overnight incubated culture. Optical den- 
sity (OD) was measured for all 96-well plates at 600 nm. Bacterial 
growth was positively evaluated when ODgoo was higher than 0.2 
(that is 5-fold higher than the untreated control) after 24 h. 

STATISTICS 

Significant differences of 16S rRNA gene based fingerprints 
between all samples were calculated with permutation analysis of 
pairwise similarities using permtest package for R statistics (The R 
Foundation for Statistical Computing Version 2.1.1) (Kropf et al., 
2004). Other significant differences were calculated using ANOVA 
with post-hoc test Scheffe in SPSS-PASW Statistics v. 18. 
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RESULTS 

MOLECULAR FINGERPRINTS OF 16S rRNA GENES OF THE BACTERIAL 
COMMUNITIES 

Molecular fingerprints were performed by SSCP analysis of 16S 
rRNA genes amplified from DNA obtained from all samples to 
gain first insight into the bacterial communities. Bacterial finger- 
prints of potting soil samples (SB-PS and PS) comprised fewer 
but more dominant bands compared to the SSCP profiles from 
WB-CS, SB-CS, and CS, where a higher number of SSCP bands 
were detected (Figure SI). According to cluster analyses, bacte- 
rial communities from SB-CS and CS significantly differed by 
78% from SB-PS and PS (P < 0.001), which clearly indicated 
the high impact of soil type and the composition of bacterial 
fingerprints. In addition, the rhizosphere effect was highly pro- 
nounced: SB-CS and SB-PS were significantly different from their 
respective bulk soils CS and PS by 72% (P < 0.029) and 52% 
(P < 0.028), respectively, (Figure 1). Differences between both 
investigated plant genotypes were also identified: the WB-CS rhi- 
zosphere was separated from all others and showed only 16% 
similarity with the SB-CS rhizosphere. Noticeably, due to different 
life cycles not only the genotype but also the plant developmental 
stage were different and could account for observed effects. 

PYROSEQUENCING-BASED 16S rRNA PROFILING OF THE BACTERIAL 
COMMUNITIES 

A pyrosequencing-based analysis of partial 16S rRNA gene 
sequences was employed to survey the diversity and composi- 
tion of the bacterial rhizosphere communities in WB-CS, SB-CS 
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FIGURE 1 | Dendrogram of 16S rRNA gene-based SSCP profiles of 
bacterial communities in association with wild beet (WB-CS) and 
sugar beet rhizosphere grown in coastal drift line (SB-CS) and 
potting soil (SB-PS) and their respective bulk soils (CS, PS) with 



and SB-PS, and the respective bulk soils CS and PS. In all sam- 
ples, we recovered between 5578 and 9567 quality sequences with 
read lengths ranging from 430 to 450 bp. Prior to further analysis, 
read numbers were normalized to 5578 for each sample. To cal- 
culate rarefaction curves and to perform taxonomic assignments, 
reads were clustered in operational taxonomic units (OTUs) at 
sequence divergences of 3% (species level). The rarefaction curves 
of the bacterial communities of WB-CS, SB-CS, SB-PS, and their 
respective bulk soils CS and PS are shown in Figure S2. At a 
dissimilarity level of 3%, the curves of WB-CS and PS sam- 
ples generally showed low slopes, but did not reach saturation. 
Accordingly, the number of observed OTUs covers only 31.4 and 
34.5% of the estimated taxonomic richness by the Chaol richness 
estimator (Table 1). That indicates evenly contributing species 
and a low number of very common or very rare species. The com- 
puted Shannon indices of diversity (H') were much higher for the 
wild beet plants grown in coastal drift line soil (WB-CS) and the 
sugar beets grown therein SB-CS, and CS (8.7, 8.1, 8.2) than for 
SB-PS and PS (6.0, 4.3). 

Altogether, 98.8% of the OTUs (total number: 5163) were 
affiliated to 10 different phyla representing at least 1% of reads 
(Figure 2): Actinobacteria and Proteobacteria were found in all 
samples. Only in WB-CS, SB-CS, and CS Acidobacteria (8.7, 
13.0, 17.1%), Chloroflexi (4.3, 7.0, 6.0%), Gemmatimonadetes 
(1.6, 4.7, 3.7%), and Planctomycetes (26.9, 39.3, 36.3%) were 
found. Bacteroidetes were found exclusively in rhizospheres 
of WB-CS (2.7%) and SB-PS (14.1%). Low percentages of 
Fibrobacteres were found exclusively in the rhizosphere of SB-CS 



four repetitions per sample. Following settings were used: dendrogram 
type: unweighted pair group method with arithmetic mean (UPGMA); 
similarity coefficient: curve based, Pearson correlation; optimization: 4%, 
position tolerance: 1%. 
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(1.0%) and Verrucomicrobia in the rhizosphere of WB-CS 
(1.7%). 

SOIL TYPE- AND PLANT GENOTYPE-SPECIFIC OTUs 

A more detailed analysis of the dominant soil and plant genotype- 
specific OTUs was performed on the base of a manual tax- 
onomic assignment by BLASTn alignment using NCBIs 16S 
rRNA gene reference database. In total, 48 OTUs were obtained 
when a cut-off level of 1% was applied (Table 2). The domi- 
nant genus within the phylum Proteobacteria was Pseudomonas 
comprising four species (OTU1726, 2094, 2281, 4102) with 
the relative abundance of 18.8% detected in the SB-PS rhi- 
zosphere, 1.8% WB-CS rhizosphere, 0.2% in the SB-CS rhi- 
zosphere, and 16.9 and 0.1% in PS and CS, respectively. The 



Table 1 | Species richness (normalized at 5578 sequences per sample) 
estimates obtained at 3% genetic dissimilarity from 454 
pyrosequencing-derived sequences of DNA extracted from 
SB-CS— Sugar beet plants cultivated in coastal drift line soil, 
SB-PS— Sugar beet plants cultivated in potting soil, WB-CS— Wild 
beet plants grown in coastal drift line soil, CS— Coastal drift line soil, 
PS— Potting bulk soil. 
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4.3 


121.3 
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8 higher number indicates more diversity; b results from the rarefaction analyses; 
c non-parametric richness estimator based on the distribution of singletons and 
doubletons. 



second major genus assigned to Stenotrophomonas comprised 
only one OTU (3731), and was detected in the SB-PS (1.8%) 
and WB-CS (0.4%) rhizosphere, and in PS (21.1%). Most of the 
Firmicutes sequences belonged to the genus Paenibacillus com- 
prising three species (OTU1033, 2251, 3065) in SB-PS (0.2%) 
and PS (8.3%), followed by Bacillus (OTU379) in the SB-PS 
(0.2%) and WB-CS rhizosphere (0.1%), and in PS (6.7%). 
Blastopirellula, the main genus within the phylum Planctomycetes, 
occurred in approximately equal abundances in the SB-CS 
(4.9%) and WB-CS (4.0%) rhizosphere, and in CS (5.2%). 
The majority of Bacteroidetes reads belonged to Flavobacterium 
(OTU254, 3902) and Pedobacter (OTU172, 3215) detected in 
the SB-PS and WB-CS rhizosphere. Within the Actinobacteria, 
Thermoleophilum occurred as the only genus (OTU745) in the 
SB-CS (3.1%) and WB-CS rhizosphere (0.6%), and CS (2.5%). 
Within the Acidobacteria the OTU945 assigned to Candidatus 
Solibacter occurred with 0.7% in SB-CS, 0.2% WB-CS, and 1.3% 
in CS. 

Comparing wild beet and sugar beet cultivated in CS, 20 
out of 38 dominant OTUs (41.7%) were detected in both rhi- 
zospheres, whereas 18 OTUs (37.5%) accounting for 2.3% of 
total reads occurred exclusively in WB-CS, i.e., OTU172, 254, 
1122, 1155, 2094, 2281, 2797, 3190, 3215, 3262, 3731, 3889, 
3902, 3915, 3972, 4395, 5140, 5143 (Table 2). The WB-CS- 
specific OTU with the highest percentage of reads (0.5%) was 
OTU2094 assigned to Pseudomonas grimontii. Remarkably, none 
of the OTUs was found to be specific for SB-CS. In the soil 
type comparison of sugar beet plants cultivated in CS and 
PS, 12 specific OTUs (25.1% of total reads) were found in 
SB-CS and 25 in SB-PS (33.7%). In SB-CS, the most dom- 
inant and specific OTU (5.6%) was assigned to Angiococcus 
disciformis (OTU1264) and in SB-PS to Enterobacter sp. with 
6.1% (OTU3915). Altogether, a high impact of soil type and 
also of plant genotype on the rhizosphere communities was 
observed. 
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FIGURE 2 | The structure of bacterial communities at phylum level of wild 
beet rhizosphere (WB-CS), coastal drift line bulk soil (CS), sugar beet plants 
cultivated in coastal drift line soil (SB-CS), sugar beet plants cultivated in 
potting soil (SB-PS), and potting bulk soil (PS). Relative composition of 



major phyla (>1 % of total reads) was determined by 454 amplicon 
pyrosequencing of 16S rRNA extracted from all samples. Multi-colored charts 
at the legend are shown for each phylum and sample correspondingly. Phyla to 
which <1 % of reads were assigned were summarized as "other." 
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Table 2 | Taxonomic classification and relative abundance of dominant OTUs (cut-off level: 1%) from rhizosphere and soil samples. 

OTU Phylum Genus Species Accession SI WB-CS CS SB-CS SB-PS PS 

945 Acidobacteria Candidatus Solibacter Candidatus Solibacter NFL074351.1 90 0.19 1.27 0.71 0.00 0.00 

usitatus Ellin6076 

745 Actinobacteria Thermoleophilum Thermoleophilum NR_036932.1 84 0.62 2.49 3.05 0.00 0.00 

minutum 

3902 Bacteroidetes Flavobacterium Flavobacterium NRJ09728.1 98 0.01 0.00 0.00 1.11 0.00 

cutihirudinis 

254 Flavobacterium NRJ085371 98 0.27 0.00 0.00 1.08 0.00 

myungsuense 

3215 Pedobacter Pedobacter borealis NR_044381.1 99 0.02 0.00 0.00 2.03 0.00 

172 Pedobacter nyackensis NR_044380.1 99 0.03 0.00 0.00 2.10 0.00 

3911 Firmicutes Alicyclobacillus Alicyclobacillus pomorum NR_024801.1 99 0.00 0.00 0.00 0.00 1.17 

379 Bacillus Bacillus ginsengihumi NR_041378.1 99 0.09 0.00 0.02 0.24 6.68 

4794 Brevibacillus Brevibacillus reuszeri NR_040982.1 99 0.00 0.00 0.00 0.02 2.28 

3224 Cohnella Cohnella formosensis NRJ09515.1 99 0.00 0.00 0.00 0.03 1.54 

1531 Oxalophagus Oxalophagus oxalicus NR_036979.1 98 0.00 0.00 0.00 0.01 1.13 

1033 Paenibacillus Paenibacillus koleovorans NR_024752.1 96 0.00 0.00 0.00 0.02 1.75 

2251 Paenibacillus sediminis NRJ08601.1 99 0.00 0.00 0.00 0.18 3.73 

3065 Paenibacillus terrigena NR_041398.1 97 0.00 0.01 0.00 0.00 2.81 

3895 Paenisporosarcina Paenisporosarcina indica NRJ08473.1 99 0.16 0.03 0.02 0.11 2.49 

2595 Tumebacillus Tumebacillus NR_043849.1 97 0.00 0.00 0.00 0.00 1.66 

permanentifrigoris 

616 Planctomycetes Blastopirellula Blastopirellula marina NR_029226.1 88 0.96 1.05 0.89 0.00 0.00 

2186 Blastopirellula marina NR_029226.1 90 1.44 2.09 1.59 0.00 0.00 

2807 Blastopirellula marina NR_029226.1 90 1.65 2.11 2.43 0.00 0.00 

2164 Proteobacteria Achromobacter Achromobacter spanius NR_025686.1 99 0.00 0.01 0.00 0.16 7.15 

5140 Aeromonas Aeromonas veronii NRJ02789.1 99 0.03 0.00 0.00 3.87 0.00 

3190 Agrobacterium Agrobacterium NR_041396.1 99 0.33 0.01 0.00 1.36 0.33 

tumefaciens 

4395 Brevundimonas Brevundimonas bullata NR_025831.1 99 0.05 0.00 0.00 1.37 0.00 

3816 Cellvibrio Cellvibrio fulvus NR_025210.1 99 0.30 0.00 0.09 1.42 0.00 

1264 Cystobacter Angiococcus disciformis NRJ04864.1 85 1.09 8.7 5.63 0.00 0.00 

2092 Angiococcus disciformis NRJ04864.1 85 ^^^^ 2.60 0.00 0.00 

4604 Desulfobacca Desulfobacca NR_074955.1 85 0.47 1.22 1.48 0.00 0.00 

acetoxidans 

3651 Desulfomonile Desulfomonile limimaris NR_025079.1 87 1.10 1.42 1.18 0.01 0.00 

3799 Desulfonema Desulfonema ishimotonii NR_025991.1 87 1.13 0.80 0.71 0.00 0.00 

4187 Desulfotalea Desulfotalea arctica NR_024949.1 84 0.29 1.10 1.30 0.00 0.00 

(Continued) 
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Table 2 | Continued 



OTU Phylum 

187 

3262 

3915 

276 

5143 

3972 
281 

3889 
3620 
2797 

2281 

1726 

2094 
4102 
1122 

1155 
3731 

3680 



Genus 

Devosia 



Duganella 



Enterobacter 



Geobacter 



Janthinobacterium 



Methylophilus 



Novosphingobium 



Paracoccus 



Pelobacter 
Phenylobacteriwn 



Pseudomonas 



Rhodopseudomonas 



Sphingopyxis 



Stenotrophomonas 



Variovorax 



Species 

Devosia crocina 



Accession SI WB-CS CS SB-CS SB-PS PS 

NR_044213.1 98 0.07 0.00 0.02 1.27 0.00 



Duganella phyllosphaerae NR_1 08529.1 



99 



0.03 



0.00 



0.00 



Enterobacter sp. 



NR_0 74777.1 



99 



0.14 



0.00 



0.00 



Geobacter grbiciae 



NR_041826. 



84 



1.08 



2.25 



2.41 



Janthinobacterium 
agaricidamnosum 



NR_026364.' 



99 



0.01 



0.00 



0.00 



Methylophilus flavus 



NRJ04519.1 



99 



0.04 



0.00 



0.00 



Novosphingobium 
resinovorum 



NR_044045.' 



99 



0.41 



0.01 



0.02 



Paracoccus limosus 



NR_1 09093.1 



99 



0.15 



0.00 



0.00 



Pelobacter acetylenicus 

Phenylobacterium 
haematophilum 



NR_029238.' 
NR_041991.1 



86 
99 



0.13 
0.03 



0.98 
0.00 



1.18 

0.00 



Pseudomonas 
brassicacearum subsp. 
brassicacearum 



NR_0 74834.1 



98 



0.12 



0.01 



0.00 



Pseudomonas 
entomophila 



NR 102854.1 



Pseudomonas grimontii NR_025102.1 



Pseudomonas poae 



NR 102514.1 



Rhodopseudomonas 
palustris 



NR_036771.1 



Sphingopyxis italica 



NR 108877.1 



Stenotrophomonas 
maltophilia 



NR 074875.1 



Variovorax paradoxus 



NR 074646.1 



1.17 



99 0.39 0.00 0.00 1.82 



0.00 



6.06 0.02 



0.00 0.00 



1.86 0.00 



1.10 0.00 



1.62 0.00 



1.00 0.00 



0.00 0.00 
1.91 0.20 



1.01 0.00 



99 1.11 0.12 0.1E 



99 0.52 0.01 0.00 2.52 0.01 



98 0.00 0.00 0.00 1.20 0.00 



98 0.03 0.00 0.00 0.23 1.85 



99 0.05 0.00 0.00 0.50 4.44 



99 0.32 0.00 0.14 1.82 0.10 



Relative abundances ( %) of OTUs in the rhizosphere and bulk soil based on pyrosequencing data. Accession numbers of the NCBI GenBank accession was indicated 
for the closest representative sequence to the OTUs with the corresponding similarity index SI in %. Relative abundances are indicated by colors ranging from low 
abundances (yellow! raising to a high abundance (red). Non-detected OTUs in respective samples are white colored. 



COMPARATIVE ANALYSES OF RHIZOSPHERE AND BULK SOIL 
RESP0NDERS 

Taxonomic comparisons were done at genus level for repre- 
sentative sequences of all OTUs obtained by clustering at 97% 
similarity. Classification was performed using the RDP naive 
Bayesian rRNA classifier and Greengenes 16S rRNA database. 
Overall, 37.4% of the sequences could be assigned to 70 gen- 
era with a confidence score of >80% (Table S2). By comparing 
CS and WB-CS, 44 taxa assigned to eight phyla were detected 
(including Other). In general, 31 taxa responded positively in the 
rhizosphere with increased ratio higher than 2, where the high- 
est ratio was found for Novosphingobium sp. (74.4) followed by 
Pseudomonas sp. (50.8). Comparing CS and SB-CS, only six taxa 
(ratio > 2) were positive rhizosphere responders with the high- 
est ratio found for Cellvibrio sp. (6.1) followed by Aquimonas sp. 



(4.3); whereas 32 taxa were not detected at all in CS. By com- 
paring PS and SB-PS seven taxa (ratio > 2) were detected with 
the highest ratio of 116.0 found in the rhizosphere and assigned 
to Sphingomonas sp. followed by Pseudomonas sp. (90.6). Fifty- 
two taxa were not detected in PS. For wild beet plants in CS and 
sugar beet plants cultivated in CS and PS, the rhizosphere effect 
could be confirmed. In general, a stronger rhizosphere effect was 
observed for sugar beet plants in potting soil compared to plants 
cultivated in coastal drift line soil. 

BACTERIAL ANTAGONISTS AND THEIR STRESS RESPONSE 

In total, 576 bacteria were isolated; 192 isolates were selected from 
WB-CS, 144 from SB-CS, and 312 from SB-PS rhizospheres. The 
lower number of isolates from WB-CS and SB-CS resulted from 
the limited amount of sample material. All bacteria were tested in 
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dual culture assays for their in vitro antagonistic activity against 
various sugar beet phytopathogens. In total, 357 bacteria were 
positive against at least one of the phytopathogens. The highest 
proportion of antagonistic bacteria originated from Kings B agar 
compared to R2A (Table 3). 

The highest percentage of antagonists were found for bac- 
teria isolated from SB-PS (R2A 15.0 ± 9.7, Kings B 57.5 ± 
11.6) against Botrytis cinerea independently from the selection 
medium. Within isolates from R2A, a higher percentage of active 
bacteria were found from SB-CS compared to WB-CS except for 
Sclerotinia sclerotiorum. Similarly, for bacteria isolated from Kings 
B, a higher percentage of active bacteria were isolated from SB- 
CS compared to WB-CS against A. alternata, S. sclerotiorum, and 
V. dahliae, except for B. cinerea and _R. solani. When cultivated in 
different soil types, isolates from SB-PS (Kings B) showed a higher 
antagonistic percentage against B. cinerea and _R. solani and less 
activity against A. alternata compared to SB-CS. 

All positively evaluated bacterial isolates (357) were further 
tested for their ability to tolerate stress including desiccation, salt, 
and reactive oxygen species caused by hydrogen peroxide and tel- 
lurite (Figure 3). In desiccation assays, the survival of the cells 
over several days was tested. In general, slightly more isolates 
obtained from the Kings B medium could be re-cultivated after 
3, 6, and 16 days when compared to R2A. After 56 days, most 
of the isolates were not able to be re-cultivated after desiccation 
(Figure 3A). In salt (NaCl) stress assays, a high percentage of iso- 
lates from WB-CS and sugar beet plants cultivated in the same soil 
(SB-CS) were able to grow in the presence of up to 6% salt. The 
majority of bacteria from all samples could deal with a maximal 
salt concentration of 6%; with the exception of SB-PS, where most 
of the isolates were unable to grow above 2% NaCl (Figure 3B). In 
the presence of reactive oxygen species (ROS) induced by hydro- 
gen peroxide and tellurite, 68.0 and 28.2%, respectively, of the 
antagonists were able to cope with the lowest concentration, and 
1.1 and 4.5%, respectively, with the highest concentration. The 
highest percentage of isolates able to deal with high ROS con- 
centrations caused by hydrogen peroxide were determined for 
isolates from WB-CS (0.1 mM: 91.7%; 0.3 mM: 53.3%; 0.5 mM: 
16.7%) selected on Kings B (Figure 3C). The maximum concen- 
tration of hydrogen peroxide at which isolates were able to grow 
was 0.7 mM. None of the isolates from SB-CS and from Kings B 



were able to grow in the presence of H2O2. Among all samples, 
only a few isolates were able to tolerate high ROS concentra- 
tions from tellurite. These isolates were obtained from SB-CS. The 
maximum concentration of tellurite at which isolates were able to 
grow was 9mg ml -1 (Figure 3D). 

Bacterial isolates with the highest antagonistic activity and 
stress tolerance were identified (Table 4). From all samples, 
isolates of Pseudomonas and Stenotrophomonas species were 
identified. Isolates from the sugar beet rhizospheres (SB- 
CS, SB-PS) were mainly identified as Staphylococcus and 
Stenotrophomonas species. In contrast, from WB rhizospheres 
more diverse bacterial species were isolated, i.e., Curtobacterium 
sp., Erwinia sp., Microbacterium sp., Micrococcus sp., Streptomyces 
sp., Sphingobacterium sp., Agrobacterium sp., Luteimonas sp., 
Pseudomonas sp., Rheinheimera sp., and Stenotrophomonas sp. 
Based on their 16S rRNA gene similarity, the abundance of bac- 
terial isolates in the amplicon libraries was calculated. In all sam- 
ples, 16S rRNA sequences were detected of species to which the 
active isolates were affiliated. Their relative abundances ranged 
from 0.01 to 2.94% with the highest abundances of Pseudomonas 
and Stenotrophomonas. 

DISCUSSION 

In this study, we found differences within the microbiome com- 
position of the sugar beet rhizosphere and its wild ancestor. 
At phylum level, the rhizosphere microbiome of wild beet and 
domesticated sugar beet plants (WB-CS, SB-CS) resembled each 
other independently from the climate and plant development, 
whereas, at a higher taxonomic resolution plant genotype-specific 
patterns were identified. This plant genotype-specific effect was 
confirmed for functional traits. The wild beet rhizosphere (WB- 
CS) was colonized by a low number of bacteria with antagonistic 
activity against pathogens but their antagonists showed a high 
potential to cope with abiotic stresses. Conversely, the sugar beet 
cultivar was able to enrich a high antagonistic potential from both 
soil types but harbored less antagonists with high stress resistance. 
In addition, an impact of both soil types (CS, PS) on the bacte- 
rial composition was found. However, the differences in structure 
and function of the microbiomes of ancient and modern beets 
underline the relationship between the plant genotypes and their 
associated bacteria and confirmed our hypotheses. 



Table 3 | Antagonistic potential toward various plant pathogens. 



Origin 


Cultivation medium 


No. of isolates 3 


A. a. (%) 


B. c. (%) 


R. s. {%) 


S. s. (%) 


V. d. (%) 


WB-CS 


R2A 


96 


7.3±10.2a 


4.2 ±7.1 a 


0.0 ±0.0 a 


2.1 ±3. 5a 


2.1 ±3.5a 


SB-CS 


R2A 


72 


12.5±3.3a 


10.4±8.9a 


2.1 ±3. 5a 


2.1 ±3. 5a 


12.5±16.8I 


SB-PS 


R2A 


120 


1.7 ±2.1 a 


15.0 ±9.7 a 


0.8±1.5a 


1.7 ±2.1 a 


n.d. 


WB-CS 


Kings B 


96 


17.7 ±15.6 a 


45.8±20.2ab 


10.4±8.9a 


1.0±1.8a 


4.2 ±2.9 a 


SB-CS 


Kings B 


72 


20.8±8.2a 


27.1 ±3. 5a 


4.2 ±7.1 a 


16.7 ± 11 .5 b 


20.8±7.1 b 


SB-PS 


Kings B 


120 


10.1 ±5. 5a 


57.5 ± 11 .6 b 


4.6±2.5a 


n.d. 


n.d. 



Legend: SB-CS — Sugar beet plants cultivated in coastal drift line soil, SB-PS — Sugar beet plants cultivated in potting soil, WB-CS — Wild beet plants grown in 
coastal drift line soil. Inhibition of typical and new emerging sugar beet phytopathogenic fungi A. a. — Alternaria alternata Nees, B. c. — Botrytis cinerea Pers., R. s. — 
Rhizoctonia solani Kuhn AG2-2IIIB, S. s. — Sclerotinia sclerotiorum Fuckel, and V. d. — Verticillium dahliae V25. n.d. — not determined. a The lower number of isolates 
from WB-CS and SB-CS resulted from the limited amount of sample material. 
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FIGURE 3 | Stress tolerance assays. Bacterial isolates able to be 
re-cultivated after desiccation (A), grown in presence of different sodium 
chloride concentrations (B), and cultivated in presence of reactive oxygen 
species caused by hydrogen peroxide (C) and tellurite (D). In total, 357 
antagonistic bacteria, previously selected from R2A or Kings B medium, were 



tested. Legend: bacteria isolated from R2A — closed symbols, and from Kings 
B — open symbols; wild beet rhizosphere (WB-CS) — circle, rhizosphere of 
sugar beet plants cultivated in coastal drift line soil (SB-CS) — square, 
rhizosphere of sugar beet plants cultivated in potting soil (SB-PS) — triangle. 
Further statistics were indicated in Tables S3A-D. 



The composition of the microbiome of both plant genotypes, 
when grown in CS, was comparable at phylum level. However, 
resolved at OTU level, based on results from molecular fin- 
gerprint analysis and amplicon pyrosequencing of 16S rRNA 
genes, plant genotype-specific OTUs were detected. The wild 
beet (WB-CS) showed the highest number of unique bands in 
bacterial fingerprints but the bacterial communities in the rhizo- 
sphere of sugar beet (SB-CS) and in bulk soil (CS) were similar 
to each other. Based on amplicon libraries, the WB-CS rhizo- 
sphere consisted of 18 unique out of 38 detected OTUs including 
Flavobacterium, Pedobacter, and Pseudomonas spp. In contrast, no 
specific OTUs could be identified for SB-CS. The low rhizosphere 
effect in SB-CS might be a result of the shorter growth period in 
addition to the different climate. The majority of the dominant 
OTUs (20) was shared by SB-CS and WB-CS and included gen- 
era of the phyla Actinobacteria ( Thermoleophilum), Proteobacteria 
{Cystobacter, Desulfomonile), and Planctomycetes (Blastopirellula). 
In all rhizospheres (CS), the myxobacterium Angiococcus disci- 
formis was found in highest relative abundances. The more diverse 
spectrum of OTUs associated with wild beet plants compared to 
SB-CS was confirmed by the elevated Shannon index, which was 
8.7 for WB-CS and 8.1 for SB-CS. The impact of the plant geno- 
type on the associated microbiota, which was clearly shown in 
our study, was dependent on the applied method and their taxo- 
nomic resolution. This genotype-specific impact was also shown 
or other crops such as maize (Peiffer et al., 2013), potato (Weinert 
et al, 2009, 201 1), and rice (Engelhard et al, 2000; Hardoim et al., 
2011). 



In addition to the plant genotype-specific community struc- 
tures, we found differences in in vitro functions of the isolates 
as well. The proportion of strains with antagonistic in vitro 
activity against phytopathogens was lower in the rhizosphere WB- 
CS in comparison to the domesticated sugar beet SB-CS and 
SB-PS. The lower antagonistic activity of wild beet-associated 
bacteria can be traced back to the fact that pathogen pressure 
barely exists in coastal drift-line soil as it is influenced by sea- 
water spray, tidal flow, and storms (Biancardi et al, 2012). In 
contrast, modern sugar beets are threatened by many diseases 
including fungal pathogens which were involved in this study 
(Khan, 2013). Under agricultural conditions, the pathogen pres- 
sure is often much higher than in natural habitats; one reason 
for this can be the naturally occurring biodiversity which pro- 
tects ecosystems from spread of diseases (Latz et al., 2012). 
B. maritima grows under extreme saltwater conditions, and 
can tolerate both high salt concentrations in soil and severe 
drought (Shaw et al., 2002). Wild beet-associated bacterial iso- 
lates showed a higher extend of stress tolerance than isolates 
of domesticated sugar beet plants with in vitro antagonistic 
activity. These isolates were identified mainly as Pseudomonas 
and Stenotrophomonas species; members of both are well-known 
for their ability to cope with biotic and abiotic stress (Alavi 
et al., 2013; Zachow et al., 2013). Within the amplicon libraries 
their relative abundances was lower than 1%. However, the pro- 
tective functions can be fulfilled by a minor fraction, which 
may contribute to important ecosystem functions (Pester et al., 
2010). 
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The community composition of both sugar beet genotypes 
showed clear soil type-specific effects. In bacterial fingerprints, 
SB-CS showed only 36% similarity to that of SB-PS. Based 
on the pyrosequencing approach, the phyla Actinobacteria and 
Planctomycetes were pin-pointed as determinants for these differ- 
ences. These taxa were found exclusively in WB-CS and SB-CS. 
Planctomycetes are known as soil oligotrophes and were enriched 
in the bulk soil compared to the maize rhizosphere (Peiffer et al., 
2013). Bacteroidetes and Proteobacteria were important members 
of the rhizosphere community. At OTU level the most domi- 
nant OTUs belong to the genus Pseudomonas. Sugar beets were 
previously reported as highly colonized by Pseudomonas species 
using a broad variety of techniques (Lambert et al., 1990; Thrane 
et al, 2000; Zachow et al, 2008; Mendes et al, 201 1). In the cur- 
rent study, Pseudomonas species were detected within all bacterial 
communities using fingerprint and pyrosequencing analysis. In 
bacterial fingerprints, Pseudomonas species were mainly detected 
in samples grown in PS and were less dominant in the rhizosphere 
samples cultivated CS. Interestingly, pseudomonads appears to 
have a minor relevance for the investigated Beta genotypes grown 
or cultivated in their natural habitat. 

In this study, we compared the microbiome of both ances- 
tral and domesticated beet rhizospheres and linked functions to 
particular isolates within the bacterial community. These results 
support the suggestion by Wissuwa et al. (2009) to apply the 
knowledge of plant genotype-specific traits of associated microor- 
ganisms for breeding strategies. Moreover, wild types of crops 
as well as soil from their original distribution area are impor- 
tant sources for isolates with antagonistic activity against plant 
pathogens or with stress protecting activity for their hosts. 
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